Vortex Collisions: Crossing or Recombination? 
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We investigate the collision of two vortex lines moving with viscous dynamics and driven towards 
each other by an applied current. Using London theory in the approach phase we observe a non- 
trivial vortex conformation producing anti-parallel segments; their attractive interaction triggers 
a violent collision. The collision region is analyzed using the time-dependent Ginzburg-Landau 
equation. While we find vortices will always recombine through exchange of segments, a crossing 
channel appears naturally through a double collision process. 

PACS numbers: 74.60.Ec, 74.60.Ge 
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The important role that vortex collisions play in the 
non-equilibrium properties of type-II superconductors 
was recognized long ago [|J0. E.g., vortices subject to 
a longitudinal (parallel to field) current are unstable to- 
wards helical expansion, leading to dissipation, but this 
instability is contained if the vortices cannot cut through 
each other. This concept goes back to Josephson [Q, 
who assumed that the flux-line cutting process was en- 
ergetically too costly to occur, a conclusion that has been 
rejected by Brandt et al. Q on the basis of the first ex- 
plicit calculation of the cutting energy. More recently, 
the extent to which entangled vortices can cut through 
each other has been argued H to determine the dynamics 
of the vortex-liquid phase, which dominates much of the 
H-T plane for high-T c superconductors [||. Assuming 
an energy barrier E/ x to vortex crossing, an exponential 
dependence r, rj cx exp(— U x /T) of the relaxation time 
t and the viscosity rj of the vortex liquid was predicted 
HQ . Both flux transformer- and c-axis transport exper- 
iments [|| are particularly suitable to test these ideas. 

In the original calculation of the crossing energy 
U x , Brandt et al. || analyzed the intersection of two 
straight vortices and calculated their configurational en- 
ergy within London and Ginzburg-Landau theory, see 
also Wagenleithner ||. Later, Sudbo and Brandt fTc[ | 
found a lower barrier when all the conformal degrees of 
freedom of an elastic line are included: the lines bend so 
as to take advantage of the attraction between oppositely 
directed vortices. The energetics of crossing in the pres- 
ence of the surrounding vortex lattice, as well as defects 
such as two-line and three-line twists, have been calcu- 
lated, both in the London model and using the lowest 
Landau level approximation E2 ]. While these previous 
studies have concentrated on metastable configurations, 
here we are concerned with the vortex scattering pro- 
cess itself and particularly with the output topology after 
the collision: vortices, being line objects, can end up in 
two alternative topological configurations, reswitched or 
crossed (see Fig. 1). In this letter we analyze the condi- 
tions selecting between these two scattering channels. 



Vortex collisions have attracted a lot of interest in a 
variety of fields, ranging from turbulence in superfluid 
Helium jl3| to galaxy formation via cosmic strings Jl4) . 
While the vortex dynamics is hamiltonian (Schrodinger 
type or massive) in these systems, the vortex dynamics 
in type II superconductors is dissipative in general, im- 
plying that a collision has to be driven by an external 
current J. Also, while the usual scattering of particles 
is characterized by a small number of parameters, the 
vortex collisions studied below involve line objects with 
an infinite number of degrees of freedom that eventually 
could influence the result of the collision. Here we study 
vortex collisions for a set of typical initial geometries and 
drives and extract the generic information at the end. 




FIG. 1. Initial configuration of two colliding vortices en- 
closing an angle 2i9o and pushed towards each other along the 
z-axis by the current density J directed parallel to the z— axis. 
The collision terminates in one of two topologically different 
configurations: (a) the vortices cut through each other; (b) 
the vortices recombine through exchange of segments. 

Below, we study the collision on two different length 
scales: at large distances, of the order of the penetration 
depth A, the vortices behave as elastic lines interacting 
through currents and a London description is valid. On 
the other hand, as the vortices approach on distances of 
order the coherence length £, the topology of the lines 
around the collision point is not well defined within the 
London scheme and we have to adopt the time-dependent 
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Ginzburg-Landau formalism jl5) (TDGL). The vortices 
then are described by a complex wave function which 
accounts for the evolution of the vortex core structure. 
Here, we consider the situation where these two scales 
are well separated, i.e., superconductors with a large 
Ginzburg-Landau parameter k = A/£ ~ 100 typical for 
high temperature superconductors. 

Our results in the London regime show that even for 
nearly parallel vortices and weak driving currents the 
two vortices always collide: the drag force acting on the 
vortex ends always over-compensates the local repulsive 
force. In addition, the induced torque twists the initially 
repulsive vortex segments towards a locally anti-parallel 
configuration with an attractive interaction, resulting in 
a violent collision. From our TDGL analysis at small 
scales we conclude that the vortices always recombine 
when the cores collide, independent of drive and of the 
angle of collision. Nevertheless, a crossing channel can 
appear in the system through a double collision mecha- 
nism: we find that for our geometry the vortex recom- 
bination naturally leads to an intermediate configuration 
which is not 'free', always enforcing a second collision. 
This second collision restores the original topology of the 
individual lines and results in a crossing configuration 
which is asymptotically free. These findings then lend 
support to an early proposal by Clem for a similar 
double-cutting scenario for colliding vortex lines. 

The initial configurations studied here are shown in 
Fig. 1: the two straight vortex lines Vi and V2 enclose 
an angle 2i9o and are driven towards each other by an 
external current density J. For large distances d > £ 
between the colliding vortices the system can be de- 
scribed within the London approximation |l7||5|]. The 
vortices are line objects parametrized by the positions 
r M (z), where fx = 1,2 for the vortices Vi and Vi- We 
concentrate on isotropic superconductors (generalization 
to the anisotropic case is done by rescaling [p.8[ ) and write 
the free energy functional in the form 



•^[{17,}] = f ^ / dr„ ■ drj 
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the first term describes the pairwise interaction between 
vortex segments, with £0 = ( < i ) o/47rA) 2 the basic energy 
scale ($0 = hc/2e is the flux unit). The terms fi = v 
correspond to the self energy of the vortices while the 
terms /i ^ v describe the interaction between them. The 
second part in (Q) accounts for the energy of the vortex 
cores; the value Cq = 0.38 is found from comparing the 
vortex energy within the London theory with the value 
derived from a Ginzburg-Landau analysis ]l9| . 

Below we concentrate on the vortex Vi; the configu- 
ration of vortex V2 follows from symmetry. The forces 



acting on each point of Vi are the Lorentz force due to 
the current density J, the friction force generated by the 
dissipation in the vortex core, and the elastic and inter- 
action forces as given by the functional derivative of the 
free energy functional ([!]) . The dimensionless equation of 
motion for a point ri on vortex Vi then takes the form 



,1=1 •'v^ 



Pu i f_M ^) ni A(dr M Ar lM ) 



PI 



c [ni A (ni A ki)] + k(J A ni), 



(2) 



= ri - r M and p M = y/\r lfl \ 2 + 1/k 2 ; ni = 
denotes the tangent vector and ki = 
2 (' = d/dz and " = d 2 /dz 2 ). Note that 



where ri„ = 
ri(z)/ri(*)| 

Eq. (E) is independent of the parameterization of the 
line; the forces acting on the segments are locally or- 
thogonal. We have chosen units: [r] = A, [Force] = £q, 
[J] = Jo = £oc/£,4>o with jo of the order of the critical 
current density, and [t] = r?;A 2 /sq with rji = 4^/2-k^ 2 c 2 p n 
the viscosity per unit length as given by Bardeen and 
Stephen ]lj| (p n denotes the normal state resistivity). 
Initially, the two vortex lines are separated by a distance 
d w A and enclose an angle 2i?o, see Fig. 1. In order to 
describe the vortex collision including all dynamical and 
configurational degrees of freedom we solve the equation 
of motion numerically using a second-order Runge-Kutta 
method. The elastic forces help stabilizing the algorithm 
provided that the distance S between neighboring points 
remains smaller than 1/k and the time step is ~ S 2 . 

The integration of the equation of motion for various 
angles $0 an d current densities J reveals two generic 
types of collisions depending on the initial angle. We 
first discuss the results for a 'steep' collision with a small 
initial angle $0 < 7r /4, see Fig. 2. As the vortices en- 
ter the range of interaction, the vortex segments close to 
the x-axis are subject to a strong repulsive interaction. 
This repulsive force is compensated by the Lorentz force 
due to the applied current density J and the drag force 
generated by the vortex ends. The latter is the result of 
the elastic forces in the line: the vortex ends far away 
from the center do not interact and are pushed by the 
Lorentz force towards the origin, their motion dragging 
the central part along towards collision. With decreas- 
ing separation the torque generated by the interaction 
grows and twists the central parts of the lines, thereby 
increasing the angle between the central segments. As 
this angle increases beyond ir/2 the interaction turns at- 
tractive, the central vortex segments turn anti-parallel 
and collide. After the collision these segments remain 
bound at a distance d < £ with both cores overlapping: 
the Lorentz force due to J is not sufficient to separate 
the anti-parallel segments until the drag provided by the 
vortex ends has become large. Note that even a small 
Lorentz force (e.g., for a small angle $0 and current drive 
J) is sufficient to trigger a collision due to the line twist. 
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In the second, 'flat' type of collision with $ > 7r /4, all 
forces drive the vortices toward each other and they col- 
lide rapidly; the collision scenario then simplifies without 
the initial repulsive phase, the central vortex segments 
twisting smoothly and colliding violently. 




FIG. 2. Vortex configurations at constant time steps for a 
'steep' collision with initial angle $o = tt/8; a driving current 
J = 0.4Jo has been chosen. Note the various behaviors of the 
central segments during the collision: repulsion, twist, attrac- 
tion, and collision with binding ((a) side view, (b) top view, 
(c) front view, units in A). 

At distances ~ £ the London model is not valid; the 
vortex cores overlap and the interaction is dominated by 
the condensation rather than electromagnetic energies. 
We then refine the discussion of the vortex collision by 
solving the time-dependent Ginzburg-Landau equations 
|L5| within a small box of size ~ 10£ around the collision 
point. The vortex configuration is described by a com- 
plex wave function ip(r, t) undergoing a time evolution 

(d t - i&ip = (V/k + iA) 2 v + (i - M 2 M (3) 

where (f> is the electrochemical potential and we use the 
units: [r] = A,JA] = V2XH C , [i/>] = y/^/fi, [t] = r 
= h"f/\a\, and [0] = h/2er. Here, a and (3 are parameters 
in the Ginzburg-Landau free-energy, H c is the thermody- 
namic critical field, and 7 is the damping coefficient, see 
|^5| . For a clean superconductor 7 = l/2fc F £o (£0 is the 
BCS coherence length) resulting in a ballistic time scale 
t = £o/ w f, while in the dirty case 7 = 3/2k F l (I the 
mean free path) with a diffusive time r = 3£g/Zu F (fc F 
and v F denote the Fermi wave vector and Fermi veloc- 
ity). Eq. (||) is completed by the Maxwell equations and 
the constitutive relations for the superconductor |15|] , of 
which we need the expressions j„ = — 2£(<9 t A + V</>/k) 
and j s = i[ip*(V/K + iA)ip — c.c] for the normal and su- 
perconducting current densities with £ = m/3/4e 2 ?ip„7; 
the continuity equation d t p + V • (j„ + j s ) = with p 
the charge density can then be used to determine the 
potential 4>. Using the values of (3 and 7 derived from 
microscopic theory Ejj] the coefficient £ takes the value 
£ = 14C(3)/tt 4 in the dirty case, while £ = 4l/ir 3 e~ c to 



for a clean superconductor (C is the Euler constant); be- 
low we concentrate on the parameter-free dirty situation. 
As we are interested in the core region (i.e., scales smaller 
than A) we may ignore transverse screening by choosing a 
gauge with A = 0. Furthermore, neglecting weak charg- 
ing effects [^5) the requirement of divergence free flow, 

£V 2 = {i/2)[ip*V 2 ip - ipV 2 4>*\, (4) 

determines the electrochemical potential <f> needed to 
push the normal current density j„ through the vortex 
core. We use the initial condition 

ip(r, t = 0) = ipiifavl — v 2 exp(— invz), (5) 

where ipi, i = 1, 2 denote the wave functions for the 
vortices Vi [ip(r,ip) — («r/v n 2 r 2 + 2) expiip] appropri- 
ately rotated and translated. The last factor in (|^) 
describes the driving current j s = 2v(l — v 2 ) along 
the z— axis with a tuning parameter v and a prefactor 
that takes into account the reduction in the superfluid 
density due to the applied current. This applied cur- 
rent is maintained during the time evolution through 
the boundary conditions d z ip — —invip on the top and 
bottom of the box. We always adapt the dimensions 
of the box such that the cores penetrate through the 
sides in order not to disturb the applied current den- 
sity; we use two types of boundary conditions on the 
box sides to make sure that our results are not influ- 
enced by our specific choice: i) the boundary condition 
Vj_V> = guarantees that no currents leave the box 
through the sides, and ii) we impose the boundary con- 
dition V>( r ? t) = 4>i ( r — v 0V'2(r + vf)v 1 — v 2 ex-p(-iKisz) 
(c.f. (j^)) with the velocity v adapted to the local veloc- 
ity of the vortex close to the boundary. We start with 
the vortices separated a distance of a few £ in order to 
have the wave functions relax before they enter the colli- 
sion region and integrate Eq. ([|) using the Forward- Euler 
method which converges well for our dissipative dynam- 
ics. The Poisson equation (|J) for $(r, t) is integrated at 
each time step using a conjugate gradient method |p0| . 

A systematic survey for different angles and applied 
current densities smaller than the depairing current pro- 
vides the following results: any collision, independent 
of initial angle or drive, makes the vortices recombinc 
through exchange of segments, see snapshots (a) to (c) 
in Fig. 3. The resulting reswitched configuration is not 
'asymptotically free': the drive current now acts differ- 
ently on the newly built vortex lines and pushes the two 
vortex ends of each line into opposite directions, inducing 
a twist of the central vortex segments. As these segments 
turn around, the Lorentz force changes direction and 
the vortices are driven towards a new collision (snapshot 
(d)). The second collision (snapshot (e)) then enforces 
another recombination, resulting in a 'crossing' topology 
as if the original colliding vortices had cut through each 
other (snapshot (f); we have found this double collision to 
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occur in all our simulations with "do/ir € [1/12, 5/12] and 
j/jo G [0.15,0.75]). We thus find a rather unexpected 
result: while each (local) collision induces a change in 
topology through recombination, double collisions occur 
naturally due to geometrical constraints and lead to vor- 
tex crossing — the asymptotic regime where the vor- 
tices can separate is only reached when the vortices have 
crossed after a second collision. 
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FIG. 3. Time evolution of a vortex collision calculated 
within TDGL theory for # = vr/4 and J = 0.6J . The time 
sequence (from left to right and top to bottom) reveals a 
double collision process: (a)-(c) first collision with switched 
segments due to recombination; (d) central segments twist 
and prepare for second collision (e) ; (f ) second recombination 
with asymptotically free vortices in the crossing topology. 

While we have studied the generic situation of two 
individual colliding vortices, vortex collisions occur and 
are important in other situations, e.g., within the vortex 
solid or liquid phase of (high-T c ) superconductors. The 
'asymptotic' regimes defining the initial and final states 
of the collision may then look quite different and the sec- 
ond collision producing vortex crossing, which is generic 
to the above geometry, may be absent. An example for 
such a situation is the expansion of a vortex twist |Tl| ] 
in a vortex lattice or liquid phase. The interaction with 



neighboring vortices inhibits the second twist of the seg- 
ments and each collision terminates after one encounter 
- as a result, the vortex twist (or ring) smoothly ex- 
pands through recombination processes. Such recombi- 
nation processes have also been suggested to relax the 
stress in driven vortex systems [Q. Furthermore, vor- 
tex recombination has been assumed to help relax the 
pitch of spiral vortices in current carrying superconduct- 
ing cylinders subjected to a longitudinal field pjj ]. 

Summarizing, colliding vortices in type-II supercon- 
ductors always recombine through exchange of segments. 
However, the line nature of the vortices can enforce a sec- 
ond collision, opening a vortex crossing channel through 
a double collision. The two cases are easily distinguished 
through careful inspection of the asymptotic constraints 
enforced on the collision. 
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